Parton picture of inelastic collisions 
at transplanckian energies 



O.V. Kancheli* 
Institute of Theoretical and Experimental Physics, 
B. Cheremushinskaya 25, 111 259 Moscow, Russia. 



Abstract 

We propose a parton model of inelastic collisions at transplanckian 
energies E S> G _1//2 , using the gravitons, whose transverse momenta are 
cut at the Planck scale, as partons. For this purpose we represent the 
gravitational shock-wave accompanying the fast particle in terms of such 
partons and take into account the higher order multiperipheral-like con- 
tributions. We argue that the internal part of this shock plane contains 
the "black" disk of radius R(E) ~ E 1 ^ 2 filled by such hard partons with 
the Planck density ~ G. When two fast particles collide the hard graviton 
production comes from the region of intersection of their black disks. The 
corresponding value of inclusive cross-section at the given rapidity and 
the impact parameter is proportional to the area of this region. The final 
state with such hard gravitons is unstable relative to the long range grav- 
itational repulsion, and this leads to the creation of the multiperipheral 
chain of black holes at later stages of the reaction. We discuss various 
details of this picture including higher order corrections and there con- 
nection to a purely classical approach; we also consider briefly possible 
changes when additional hidden dimensions are present. 



1 Introduction 



The dynamics of particles interaction at transplanckian energies is rather poorly 
understood up to now, even qualitatively. It is usually believed that at such 
energies the dominant mechanism of the interaction is a gravitational one, so 
that its contribution highly overcomes that from other interactions. There is a 
lot of papers where this problem is considered from a different direction - here 
we mention only some of them [Q - || . The renewal of interest to this topics 
in last years - is probably connected with the discovery of various brane 
world like scenarios with TeV-scale gravity jL7|. If this will be realized, then 
gravitation mediated effects can be found already at future accelerator energies. 

When two particles with high center-of-mass system (c.m.s.) energy y/s 3> 
G -1 / 2 = m p and impact parameter B < ^/s ra~ 2 come close one to another, 
then the high energy (mass) = \/s is concentrated for a small time in the 
compact region. For this energy the 'static' radius of the horizon of a black 
hole (BH) with the mass y/s is much larger than the size of the region, where 
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the energy is concentrated. Then, naively, one can expect that in the process 
of the collision (or even before), an event horizon "surrounding" this region will 
form, and after this horizon will gradually (at time ~ y^s/m^) transform to the 
boundary of a real BH with the mass ~ ^/s. But, from the other hand, the 
colliding particles arc highly ultrarelativistic and it is unclear if it is enough 
time that such a horizon can appear in a causal way. 

To refine (or reject) these conclusions number of approaches was developed. 
The gravitational field of a fast particle with energy E is concentrated in a shock- 
wave plane of thickness ~ 1 / E and is usually represented by the Aihelburg-Sexel 
(AS) metric fL8| , The process of a classical collision of two such AS waves is 
considered in the number of papers. In details it was firstly done in [[l] where 
the transition of one AS wave through another is considered in the geometrical 
optics approximation. It follows that due to the curving and later focusing of AS 
planes the curvature singularities can appear. They signal the possible creation 
of a BH in process. It is also mentioned in B,S that a numerical solution of 
the Einstein equation for this process also point to a possible appearance of 
curvature singularities. The other (but in principle the close) method jl5| is 
based on a construction of a maximally large trapped surface in the AS disks 
collision process and perhaps confirms the same conclusions. 

Unfortunately the transition of the one AS disk through another is a very 
complicated nonlinear process with the unstable behaviour, and it is unclear if 
one can believe in answers based on the geometrical optics approximation. The 
other trouble of the simple classical approach is that the curvature in essential 
parts of AS-disks is too high (3> ml) - and thus all higher order corrections 
in curvature to the effective Einstein Lagrangian can be of the same order, and 
this, like in the case of the string interaction, can change the situation drastically 
- in particular it can soften high curvature effects and freeze them on the ~ m p 
scale. Then all the AS-disk focusing picture [|l] can change for B < y/s/m,p, 
because the different degrees of freedom are relevant. 

The different approach to transplanckian collisions is the perturbation the- 
ory. The diagram with a one graviton exchange leads to the cross-section fast 
growing with the energy 

MB < B ) = J £ dk\ ~ GVfc- m ~ B^f , (1) 

where k\ min = B^ 1 is the infrared cutoff. The expression (|l|) hardly violates 
the unitarity for s/rrip 3> 1 because the allowed maximum of <r e i(B < Bq) ~ 
Bq. The eikonalization of the one graviton exchange cures it, and leads || to 
the elastic amplitude 

A(s, b) = 1 - e lS , 5 = 8Gs In ^ ■ (2) 

b 

This amplitude has a simple physical interpretation at the lab. frame of the one 
of colliding particles, where the phase <5(s, b) coincides with the time delay of 
the target ("test") particle when it penetrates through the AS shock front at 
the impact parameter b. 

Other high order corrections to (||) are also considered in the literature - they 
include the multiperipheral diagrams with gravitons and some generalizations, 
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similar to that formulated in the gluon (BFKL) case j|, ||. But no conclusions 
going far beyond that given by ([j]) and (Q) is reached so far this way. At the 
same time it is difficult to apply directly such perturbative methods to the region 
of "small" impact parameters b < mp 1 (s/m, p ) 1 ^ 4 , where, as one can expect, 
the most interesting effects take place. One needs sertan indirect (bypass) way 
which at the same time preserves the main aspects of the high order perturbation 
theory, encoded in multiperipheral diagrams and reggeons. 

In this paper we speculate how the transplanckian particle collision process 
can look out in the parton picture. One can hope that by this approach one can 
take into account the perturbative graviton mechanisms and also include implic- 
itly essential nonperturbative and string effects. The parton picture represents 
nicely all the main aspects of high energy hadronic interactions, it also shows 
and "explains" the space-time picture of the process. It can naturally take into 
account all the regge specifics of the interaction and also that one contained in 
the dual models. Therefore one can expect that similar methods can be helpful 
for gravity induced high energy interactions, partially because there is no big 
difference between graviton and gluon degrees of freedom on the string scale. 
If partons are chosen successfully, this can give a full qualitative picture of the 
process and help to find a more accurate way to consider such a phenomenon. 

Our "main partons here are the gravitons with the mean virtuality ~ m p . 
We suppose that gravitons with higher virtualities are suppressed in a wave 
function of the fast particle by some string-like mechanism Fl. And partons with 
smaller virtualities interact much slowly and lead only to corrections and to 
infrared effects. 

The outline of paper is as follows. 

Section 2. We consider some classical aspects of the collision relevant to 
our consideration. We also discuss the popular opinion that the transplanck- 
ian gravitational interaction at not large impact parameters can be adequately 
described in classical terms . 

Section 3. Here we consider how the wave function of the fast particle can 
look in terms of such partons-gravitons. We start from the Weizsacker- Williams 
(WW) like partons, decomposing the classical AS field and then try to refine 
the model by the inclusion of graviton interactions and cascading. The hard 
parton-graviton density is large in a transverse disk with the radius ~ m~ 
and we suppose that inside of such a disk the 2D parton density is saturated at 
the Planck scale. 

Section 4. We consider the collision of two such disks and study main inelastic 
processes. The resulting final state depends essentially on the impact parameter 
of the collision and consists from layers multiperipheraly distributed in rapidity 
and filled with gravitons. In every such layer the relative velocities of gravitons 
are locally thcrmalizcd, with relative energies ~ m p . The process of the creation 
of this final state takes a long time ~ m~ 2 i/s, and the order in which these layers 
are created is system dependent, like in the usual multiperipheral picture. 
Section 5. Then we consider the future time evolution of such a multiparticle 
state, created in a parton disks collision. The classical instabilities, due to a 
long range attraction between particles with not large relative velocities, can 
take place and lead to the creation of the "multiperipheral" chain of BH from 

1 The string-bit model jscj, |2l| is in some respects similar to the hard parton constructions 
we used above, but it starts directly from a fast string state 
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particles collapsing in every layer. 

Section 6. In conclusion we discuss some questions related to this picture. 



2 Some classical aspects of collision 

The gravitational field of a fast particle with the mass m <C m p = G~ 1//2 can be 
found by boosting of the static Newton field. For a very high energy E 3> m p 
it gives approximately the same as the boosted Schwarzschild field 



(0) _ IGP^Pv ^ _ , 

yfx\ +7 2 (z - (3t) 2 



where gffl is the Minkowski metric, - fast particle 4-momenta, 7 = E/m, 
(3 = p/E . By going to the limit /3 — ► 1 and making an additional singular 
in t, z coordinate transformation, one usually represents this metric in the AS 
form 

ds 2 = dx + dx~ + AGE\nx\5{x-)(dx-f - {dx±_f , (4) 

where x + = v = t + z, x~ = u = t — z, 11 are the light-cone coordinates. 
Another method to become the AS metric which can be simply applied also to 
a distributed light-like matter (~ partons), is based on the old observation that 
for a general right-moved plane metric 

ds 2 = dx + dx~ + f(x~,x±) (dx~) 2 — (dx±) 2 (5) 

Einstein equations reduce to a simple Poisson form 

R— = dif(x-, Xl _) = GT__(ar,:r ± ) (6) 

For a bunch of massless particles with energies e„ and coordinates (z n , x n ±_), 
moving along the z axis, with the energy-momentum tensor 

T — (u,xj_) = ^e n 5 2 (x±- x n ±_)5(u- Zi) . (7) 

71 

it follows, due to linearity of (0), that 

/ = 4GEY, (|r)<5(z~ ~ Zn)H*± - *n±) 2 ■ (8) 

n 

For n — 1 this reduces to the AS metric (Q). For one point-like particle the 
curvature fields are concentrated in a plane shock-front 5{x~), like the electric 
and the magnetic fields for a fast Coulomb particle, where the same structure 
5(x~) lnx 2 ± enters the expression for potentials []. For the metric of type (||) 
the only nonzero components of the curvature tensor are 

R-±-± = - d±d±f(x-,x±_) - 5(x~)xj_x±/x 4 L , (9) 

where the last expression is for a point-like particle. The curvature tensor for 
(||) has approximately the same structure as (^), but is evidently smeared in x~ 

2 The often mentioned analogy of the AS type metric with shock waves is very superficial. 
In usual shock waves the energy is concentrated not in the front discontinuity but in a long 
layer of the gas moving behind the front. 
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on mx±/E and, what is in a certain sense essential, its Weyl part belongs to a 
different canonical type. 

In connection with the previous is also such a property of the AS metric 
(|]), and also of the more general one (||), that the corresponding symmetrical 
energy-momentum pseudotensor Q^J^ vanishes - so that it is hard to define in 
a consistent way the energy contained in AS type field ^. But for the boosted 
Schwarzschild metric (^|) the is nonzero - it has only one large component 



JLL) 



G~f 2 m 2 



(10) 



as is normal for an ultrarelativistic object. The energy flow in this case is 

oo 

jll) 1 f m\f E 



e(x±) 



dx~e\ 



+ (-)(" 



and the quantity 



oo 

I d 2 xj_e(xj_) ~ (—) E 
J \m p J 



l/m 

gives the energy in the coherent gravitational field of fast particles. 

The massless test particles colliding with the AS disk on the distance x± from 
the center will be captured (trapped) by a disk for time = 8GE ln(L/x±), and 
during this time transported backward in z to the distance 8GEhiL/x± and 
then released from opposite side of the AS disk (See Fig. 1). After that particles 




Fig.l : Trajectories of "test" particles a\ and ai colliding with 
the AS disk accompanying a fast transplanckian particle. On the 
sections C\ and C2 the particles are "transported" by the AS disk. 

continue there motion with the same momenta as before the collision 



'This takes place not only for a Landau-Lifshitz form of 



n(LL) 



The other popular forms 



of for example the canonical one, also give zero. This zero is the result of the change of 
the canonical Weyl tensor type in the AS case - which by itself follows after the singular 
coordinate transformation made while going from (j^) to (|i ) . 

4 Such a behavior is a simple reflection of that phenomena that in the gravitational field of 
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Trajectories of these test particles look discontinuous (if we "forget about 
sections C\ and C2 on Fig.l, where particles simply move backward glued to the 
AS plane). One can perform the specific coordinate transformation [Q, singular 
in the AS plane, so to make the test particles trajectories, crossing the AS 
plane, continuous. Making it one cuts out a region of the space-time around the 
AS plane (of 1 /E thickness) , after that shifting and gluing together the distant 
planes. But this procedure looks very strange and risky for a nonsingular metric 
, for which the test particle trajectories are continuous. 
The classical collision of two AS disks is considered in the number of works 
but no exact solution showing how they move one through another was found. 
The most detailed analysis of this process was first given in the paper of D'Eeath 
[0 (See also ||] and papers cited there). One can mention three main approaches 
to this problem. One is based on a numerical solution of Einstein equation, 
another on the geometrical optics approximation, and the third on the analysis 
of the possible appearance of trapped surfaces during the process of the two AS 
disks collision. 

In the geometrical optics approximation one of the AS surface is treated as 
a system of many infinitesimal test particles (directrix) which cross the another 
AD disk and do not disturb them during the traverse period. The same is 
symmetrically done with another AS-disk. Because the time-delay depends on 
x±, during such a transition of the AS fronts one through another, they curve. 
At later times these curved AS disks focus on caustic surfaces, where their 
curvature becomes very big, and this signals about the creation of some BH like 
objects, possibly surrounded by horizons. As mentioned in the consideration 
of the same problem by means of a numerical solution of Einstein equations gives 
the reminding picture. 

Another approach is based on a search of possible trapped surfaces shortly 
before the AS-disk's collide (I], [l5| |l6|. If such a surfaces exist, one usually 
predicts that a black holes will be created in process of the collision and even 
estimates their mass and the corresponding cross-section for the inelastic interac- 
tion. Such a consideration is usually done in specially transformed coordinates, 
singular with respect to the "normal" one and chosen in such a way to make the 
test particles trajectories "continuous", when they cross the AS surface. This 
is a risky procedure because some effects can arise due to the singular structure 
of AS planes and the real answer can depend on that what takes place when 
two singular AS planes go one through another in the period of their collision, 
where we in fact do not know the metric. 

At the same time one can try to consider this in "usual" nontransformed 
coordinates. The picture is illustrated by Fig. 2 where the configuration of two 
disks is shown for some time before the collision. The region T is defined in 
such a way that the massless test particle can not escape from this region before 
the AS disks coincide. The space between disks is flat. But when a particle 
will try to escape from the region T it will collide with one of the disks and 
will be absorbed (glued) and transported with it towards another AS disk. The 

a "standing" particle time is delayed. So moving in this field massless particles are slightly 
slow down (remain behind in z) near the mass on time St ~ AGmln(L/xj_), where L 2> x± 
is the full distance in 2 that particles traveled. As a result they remain behind in z from 
the similar trajectories with larger Xj_ If we boost this time delay by gamma-factor ~ E/m 
we come to the space-time delay AGEln(L/x±), and this corresponds to the picture of the 
motion of the test particle captured by the AS front. 
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Fig. 2 : Particles from T - the filled region between approaching 
one to another AS disks can not escape to infinity. They collide 
with disks AS± or AS2 and after that are transported to the plane 
z = where two disks coincide(merge), before particle come out 
from other side of disks. Particles a\ and a 2 shown on Fig. are 
trapped and particle 0,3 have time to cross the ASi disk and es- 
cape. This type condition in fact fully defines the shape of region 
r for all B and y. 

longitudinal and the transverse size of T can be defined from the condition that 
the particle glued to one of ASi disks will not be released (from another side 
of the disk) before disks collide. The region T is maximal in c.m.s. and for 
B = 0, when it is cylinder with the radius x± and the height Z = 2x± = 8GE. 
In the frame, where particles have energies E±, E2 but again B = 0, the region 
r is a truncated cone with foundations Ri = 4GEi and Z = 4G(Ei + E2). For 
B 7^ the shape of T is deformed (sec p6[). 

But how can one conclude only from the existence of a region T with such a 
properties that anything singular happens in the future evolution of the AS disk 
system? Evidently that this future depends only on that what takes place on a 
complicated nonlinear stage when two AS disks (smeared in z by ~ V-Ei) coin- 
cide. And the existence of the empty region T nothing changes in this respect. 
When regarded from such a position, the trapped surface found in coordinates 
with continuous trajectories can be the artifact of using the singular coordi- 
nates. Another risky moment of using the trapped surfaces for the prediction 
of the future system evolution is that these trapped surfaces are isolated from 
the "rest world" by the the AS shock fronts where the curvature is singular (or 
partially coincide with the ASi planes). And the last moment. If for colliding 
particles we use the regular metric ([}[), instead of the AS, then the previous con- 
sideration of r (like in Fig. 2) remains unchanged. But if some trapped surfaces 
exist in this case is rather unclear now and needs an additional investigation. 
It is possible that the classical "singular" methods as used in [Q, ||, [l(| lead to 
more or less adequate answer for the BH production, but the physics of this is 
not evident. 

At the end of this section two additional remarks 

• The whole picture of the AS disks collision should be longitudinally boost 
covariant, when we change the energies of colliding particles E\ — > Ei£, 
E2 — > E2/S, leaving s ~ 1E\E2 and the impact parameter B invariant. 
By this transformation we can always select a frame in such a way that 
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it is close to the "laboratory" frame, when one of the energies is very 
big Ei ^S> nip and the another is "arbitrary" small Eq, <§; m p . Moreover, 
because the mass of the particles is "in our hands" we can make the energy 
Ei so small that it can be considered as a "test particle" with the respect 
to the particle E\. In such a system the "trapping" region T disappears 
and all the space-time before the ASi disk corresponds to the region of 
the intersection of two disks in c.m.s. (z = in Fig. 2). Now we have only 
one AS disk and an additional light-like test particles can move to infinity 
without any intersection with the ASi plane. All this means that if we 
hope that some classical curvature singularities must necessary appear in 
such a collision (remember that s 3> Trip), then, in this system, they must 
originate from an instability of the individual AS disk with the respect to 
a small disturbances. 

• The curvature in the AS shock front is high - for the transplanckian case it 
is much larger than the critical one - for some components of the curvature 
tensor l-R^Ao-l 2> Trip. Therefore in the effective gravitational Lagrangian 
the higher order terms in R^xa can be of the same order of the magnitude 
(or even bigger) than the standard R term. Such terms can be induced by 
the quantum fluctuations or "come from strings" . For the metric of type 
(|j) this will correspond to additional terms at the left hand side of Eq. (^) : 

dif(x-,x ± ) + c 2 m p 2 (dlf(x-,x ± )) 2 + ... = GT__ (11) 

These terms can radically change the structure of AS solution for E 3> m p , 
smoothing the entering (^) singulariti es p| Also, depending from signs and 
relative values of nonlinear terms in (|1 1|) , this can lead, for example, to a 
saturation of values of / in the internal parts of AS disk, where / in the 
linear case can be arbitrary big. 



3 Parton structure of fast particle with E m p 

To understand the AS shock fronts interaction, when the classical picture is 
unreliable, we introduce the "microscopic" model of the AS disk in terms of 
parton degrees of freedom. There is a number of possibilities for a choice of 
constituents (partons) Various degrees of freedom or there combinations can 
be used as partons but it is very difficult to choice between them the most 
adequate one. It is possible the string coordinates (or string bit's p0|) are the 
most appropriate ones. But the dynamics at the high string density is also 
unclear - it is possible here again we come to a particle-like media. We use the 
transverse gravitons as partons, and regulate their spectra at k± > m p in a way 
suggested by "strings" . Going to Fourier components of the metric 

]Ta A (fc)e^ = w f d 3 x e ik * ( g^{x) -g$(x) ) , (12) 

A J 

where e$ -gravitons polarization tensors, uj = + k\ , and substituting 
into the right hand side of (|lj) the AS metric components from (^) or from (||) 

5 For one light particle state one can always, by moving to a rest frame, make the contri- 
bution of additional terms in (hlh small. But for colliding transplanck particles these terms 
can be very essential, especially in a AS-disk's intersection region. 
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we find: 



From the transversality condition k^e^J — in the gauge = we 

have the following relation between longitudinal and transverse projections of 
polarization tensor : 

& - «^ ■ ("J 

With the fields normalization used in (jl^) , the full parton (graviton) number in 
a fast particle field is given by : 



N = J d 3 kn(k) = J^J2 « A ( fc ) a+A ( fc ) • (15) 



Then, substituting ( |T3| ) and (|T^) in ([l5]), we can extract the transverse graviton 
density in the AS disk : 

(a x ) 2 

dn- L (E,u,k±) ~ — '— du;d 2 k 1 _ ~ 

- (-) 2 ^« fc i - (-) a T**x , (16) 

\m p / (fe ) \m p J oj } 

These spectra differs crucially from the QED vector parton spectra by the w 
and k± dependence and by the overall growth of the parton density ~ E 2 . 

In a "full" theory one can expect that the spectra (|l6|) are essentially sup- 
pressed at k± > m p due to various nonlocal 'string' effects. Here we simply 
suppose that they are cut at fcj_ > K ~ m p . Than the following integrated over 
fcj_ spectra 

min(u>,«) 

f / E \ 2 k 2 
dn^cj) = dw- / d 2 fe i .ji- L (w,fcx) ~ ( — J ^77-2 2\ (17) 







m p ) oj(lj 2 



can be approximately used for all u). 

Now let us discuss how essential are the higher order corrections to these 
primary parton spectra. The main higher order diagrams at high energies, 
which "follow" the WW contribution - are the multiperipheral ones. In terms 
of partons they correspond to a parton cascading which increases with E the full 
number of low energy partons, and as a result the cross-sections []. Such graviton 
ladder corrections for a scattering amplitude are considered in a number of 
papers || in the pure gravity, supergravity and also for "multiperipheral" string 
amplitudes Q. We will not review here these calculation - for our qualitative 
consideration we only mention some aspects essential for us. Firstly we suppose 
that such ladders are supplemented by some vertex factors T(ki±/m 2 ) that 

6 The parton spectra (|l6[) can be found also from the WW like factorization of one graviton 
exchange diagram as is usually done for photons in QED books. It is instructive to compare 
the expression (|16j) with the general form of WW spectra for the spin J field 

E\ 2J /uiduj\ d 2 k ± 2 dx d 2 k i _ 



(k±EY J fujduj\ d z k± „ 



fc 2(2-J) 
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B k L 



Fig. 3 : The fast particle A with its parton cloud interacts with 
the target B ; (a) and (b) represent the WW amplitude and one 
cell ladder corresponding to a WW cross-section ; diagrams (c) 
and (d) represent the graviton cascading and the ladder amplitude 
corresponding to a cross-section of the interaction with the target 



smoothly cuts high virtual transverse momenta > m p on exchanged lines 
(see Fig. 3). Then for integrations over kj± in all cells of the ladder the mean 
contribution comes from the region near this upper bound fc,-j_ ~ m p . 

The averaged picture of the "ladder" -cascading process in parton terms can 
be represented as a chain of L step convolutions of the "primary" WW spectra 
n ± (u}i,Wi + i,k±) given by (|l6|). This gives: 

J duj\dPk\ ^(E^ivijki) J did2d 2 k 2 n x (wi, uj 2 , k 2 ) ... 

... J diOL-id 2 kL-i ^(wn, lul, kL) ~ (18) 

A L (E_ \ 2 J_ l-i(E_\ 
{L-l)\ \m p ) u\ \uj 1 

where 

A = „-/Ar*(|), 

and where, as have mentioned above, we included the effective graviton emission 
vertexes T in the density n 1 - to take into account nonlocal effects responsible for 
the cutoff of high kj_ . Summing over the number of cascade steps L we come to 
the "reggeized" graviton-parton spectra 

dnHE^k ± ) = (—) 2 (-T-,d 2 ^ . (19) 

where gj is the corresponding coupling constant : g± <-> e qe( i, 32 *-* vG and x - Feynman 
scaling fraction uj/E. 

7 The parton cascading gives the shift of the reggeon intercept a(0) from their bare values. 
This corresponds to an additional growth of the parton density by the factor ~ (_E/oj) a , 
when we go to smaller parton energies u). The average interval in rapidity between the parton 
splitting (ladder rungs) is ~ 1/A. In the BFKL [E3| case (vector partons) the shift of the 
pomeron intercept ct(0) from their bare value is A ~ 9qcd' ' 3u * gl uon intercept probably 
does not shifts because of the gauge nature of gluon. The gravitons intercept probably also 
do not change, by_the same reason. But the intercept of the bare 2G trajectory ( =3) 
corresponding to (|lq) can shift, if nothing prevents. 
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Such corrected parton spectra are even more concentrated £| at low ui than 
the primary WW spectra (plf), and the value of A gives shift of bare WW 
intercept ( =3). 

The other important higher order corrections, besides the parton cascad- 
ing and various rescatterings (like eikonalization) , correspond to the parton 
recombination which become essential when the parton density is high, and this 
process can lead to the parton density saturation. 

Due to the masslessnes of gravitons there is also the 'infrared' contribution, 
coming from small — > 0. In the first WW approximation this corresponds 
to the large impact parameter elastic scattering. If some cells, in hard ladder, 
contain small ki± — > than in the impact parameter space this will correspond 
to big 8b ~ k~2 steps, "transporting" part of the hard ladder far in b. In terms 
of t-channel regge amplitudes this corresponds to the k\ In k\ type singularities, 
coming from such infrared cells. 

Firstly we illustrate this in terms of the parton cascading. Suppose that the 
first step (cell) is infrared and consider it. Then the parton distribution ([l6|) at 
large transverse distances x± 3> m" 1 can be represented in the form 



dn ( E \ 2 1 - 



duid 2 x±d 2 x± \m p J u> 



,3 



where 5 is narrow ^-like function, having width of peak ~ 1. Integrating ( p0[ ) 
over k± we become the distribution of these partons in the transverse (impact 
parameter) plane x± 

d 3 n ( E \2 1 1 



dud 2 



x± 



( E \2 1 1 

— — — . 21 

V m„ / uj j x \ 



which we use for x± 3> nil . All this corresponds simply to the approximation 
dn/dxj_ ~ dn/dkj 1 . But if we insert two or more infrared cells to the cascade 
chain (|I§|), we come to a much smaller contribution at large b. Indeed, in this 
case we have the convolution of two distributions (]2l]) containing integrals 

d 2 x ± /(b 1 -x ± ) i (x ± -b 2 ) i ~ m p /(& 1 -6 2 ) 4 

\x±—bi | >m~ 1 

which for (pi — b2) 2 S> rn~ 2 gets the main contribution near the integration 
ends, where one of these cells enters in the hard regime. The distribution ( pl[ ) 
at x± 3> m" 1 refers to soft gravitons, but inserted in the cascading chain 
these partons are the sources for the next hard chain sections. Then the simple 
generalization of ( [i"8| ) and ( pl| ) gives the x± distribution 

d 3 n /E\ 2 /E\ A 1 1 dn / E \2+a m 2 

dujd 2 x± \m p J \ojJ oj 3 x\ dx\ \m p J (m p x±) i 



8 The value of A depends strongly on the behaviour of the nonlocal factor T. In the 
perturbative string approach the value of A can be small ~ , where g 3 is the string constant, 
if the "cutoff" in kj_ is on the string scale and not on the Planck scale, and the higher string 
modes entering loops will not compensate this. But if we include in cascading all other 
(not only ladder) diagrams, the value of A can be turn out to zero. This follows from the 
momentum sum rules for the parton distribution n 1 - . The additional arguments in favor of 
the condition for a final A = are given in the end of section 4 
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for hard partons with < k\ >~ nip and ii >,> rn^ 1 . 

Now consider the interaction of a fast particle, with parton spectrum (|l9|), 
with a pointlike target. The cross-section is 

2+A 



dud i k±n(E,u,k±) a(u,kx) ~ s J+ , (23) 

where a is the parton cross-section on a local target. Such a behaviour corre- 
sponds to a regge pole with the positive signature in a vacuum channel and with 
the intercept 

a(0) = 3 + A . (24) 

This regge pole can be considered as an analog of the pomeron for the gravita- 
tional interaction p| (later, for abbreviation, we will call this reggeon 2G). The 
same cross-section (|23|), written in the impact parameter space, looks like 

(^) 2+4 _J_ (26 ) 

The singularities of ( p5|) in t = —q\ can be represented by the following expres- 
sion 

e- Q '^ lns (a + a iq \\n\ + a 2 (g 2 ln^) 2 + ... ) , (26) 

giving the main factors for an associated t-channel amplitude, where terms with 
coefficients a, correspond to a contribution with i infrared cells (loops) in the 
parton cascade. These terms lead to the 6 3> m^ 1 behavior generalizing (p5|): 



m ^ ' ^ ( nip ) ( In s ^ ( 4a' In s 
a\ a2 Info 

+ Jbn^ + Jfon^f 4 



(27) 



We see that at large b only terms with one infrared cell are essential. And 
thus terms in (|26|), singular in g^, can be included in vertices connected with 
the exchange by the 2G reggeon, with the intercept (0) . On the contrary in 
massive theories like the QCD, only a first term in (|27|) is present at all b (this 
is a standard regge pole contribution). As a result this leads to a fast cutoff of 
Uin(s, b) in b and, after the unitarization, to the Froissart-like behavior ~ In 2 s 
of cross-sections. 

The cross-sections <Ti n (s,b) grow fast with E and for some E exceed the 
maximal value = 1 allowed by unitarity, when all incoming particles flux at 
given b is absorbed. This means that the parton screening becomes essential. 
The simplest way to take screening corrections into account is to use the eikon- 
alization : 

a in (s,b) = l-|S(s,6)| 2 = 1 - exp(-2aO„( S ,6) ) , (28) 



9 It seems that there is no such state between the perturbative states in known sectors of 
string theories containing the gravity. But it can be a composite state seen only at the strong 
coupling. It is possible as well that this state is masked by the singularity corresponding 
directly to a black disk (some brane-like object ?). 
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where af n is given by (E7]). The S- matrix, entering (|2cf), contains, among others, 
a "black disk" , where the parton density exceeds the "critical" one ~ m~ 2 . This 
is in the internal part of the AS disk where 

E \ !+ A /2 



< Rl(E) = mf (— ) (29) 



Outside this 'black' part where x± > R±(E) the density decreases as : 



dn 2 fR ± (E)\ 



4 



For 6 < Rj_(E) the screening of partons in a process of the interaction is es- 
sential and the expression takes it into account. But, at the same time, it 
corresponds to that all these partons are nevertheless present in the incoming 
state. And therefore, there 2 D density (in parton state) anyway highly exceeds 
the Planck density m~ 2 . But then the mean transverse momenta of these par- 
tons, due to their strong interactions, also can be much greater than the Planck 
scale. And this is inconsistent with our previous assumption that all much 
higher than m p , are cut. Therefore we also suppose that the parton density 
saturates at the Planck scale and does not increase with E in an internal parts 
of the AS disk, where xj_ < Rj_(E). The dynamical mechanism responsible 
for such a density stabilization (saturation) can be probably represented as a 
parton recombination, like that in the vector (BFKL) case 0. 
The full 5-matrix corresponding to ( pcf ) 

S(y,b) = e , Sfab) ~ a°„( s ,6) + ••• (31) 

contains also a real part of the phase Sn(y, b), which for very large b coincides 
with the "classical" AS phase 

B 



SrM „ ( ■ ) „4 



and is responsible for a large b elastic scattering. This S-matrix gives the hard 
inelastic cross-section 



^(1-e-MM) ~ nR l {y) „ m; 2 (-^) , (32) 



corresponding to the same structure of a black absorbing disk with the radius 
R± (y) ~ \fs I m p as given by the expression (^9|) . Such a disk remains a Froissart 
black disk and the method by which it is often introduced from the cikonalization 

10 There are signs that the string system undergoes some pha se transition at Hagedorn 
temperatures to some new phase, and there are arguments |19| that this phase is not a 
string like. For the high string densities the situation is probaEly the same. And more, 
there exists a popular opinion formulated in various forms that virtual (parton) string states 
with the transplanckian density of degrees of freedom 3> 'i~rip D (even local) are inaccessible, 
because their contribution is strongly damped (possibly exponentially in density) in the wave 
function. Despite the fact that such a reduction of high density states can be encoded in very 
general constructions, like holographic principle or M-theory, the concrete mechanism doing 
this in all cases can be simply a more fast recombination of additional degrees of freedom than 
their creation, when some critical density is reached. 
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of supercritical pomeron exchange. The difference is the behaviour of the black 
disk radius with the energy : R ~ mT 1 ln(s/m 2 ) for the Froissart case (for 
massive theories), and R ~ m~ 2 v / s for the massless gravitation 

We considered above only partons with u ~ m p , where the majority of hard 
partons are concentrated. We can extend this to partons with different lj and 
find their distributions. Simple generalization of it gives the parton spectra at 
arbitrary lj and x± 

dn = h(E,u>,x±) (mld 2 x±) — , h = (— ") -, — -j , (33) 

u> \lj/ (xj_m p ) 

where the density n is dimensionless and boost invariant. This distribution takes 
an even more simple form if we go to rapidities Y — lnE/m p , y ~ \yilj /m p and 
to £ = ln(a;j^m p ) 2 : 

dn = e (2+A)(ir-«)-c dy dC (34) 
From the condition fi ~ 1 wc find the black disk radius for various u> : 

R 2 ± (E,u>) = m7 p 2 (-J (35) 

Thus the saturated part of a hard parton configuration in (£, y) space consists 
from a "cone", with the boundary £ < (2 + A)(Y — y) (this corresponds to 
x±_ < R^EjUi)), filled with the density ~ 1 (~ m p 3 in units x±). This 
cone is surrounded by a more "diffuse" soft (lj -C m p ) cloud, where the parton 
distribution is given directly by (|34|). 

For lj <C ni p the parton density n in this soft tail can be very large n^S> 1. The 
mean transverse localization of these partons is Sx± ~ 1/lj. Comparing it with 
the size of the region R±(E,lj) where h > 1 we can find the soft critical radius 
and the parton frequency 

/ E \ l+2A/(2-A) 

R(E)~m-\—) , Q~1/R(E), 

which define the border of the fully coherent parton region. For x± > R(E) the 
parton interaction with the target is week and is represented by the "classical" 
elastic scattering and the soft classical gravitational wave emission. 



4 Collision of two AS disks 

Firstly, let us consider an interaction of two transplanckian particles, "repre- 
sented" by such parton disks, in a formal way - by the regge like machinery. 

11 There is an interesting question, concerning the transparency of this black disk, filled with 
hard gravitons at a saturated density, and of the dependence of this transparency from energy. 
At first it seems that for the saturated parton disk, with the finite 2-dimensional density of 
degrees of freedom, the transparency (= \S(y,b)\ 2 ) should remain finite (grey disk) and not 
— > when E — > oo. But this contradicts probably to the longitudinal boost invariance of the 
transparency at a given impact parameter calculated in the parton model. The same problem 
emerges in the QCD when we try to estimate the transparency of the saturated Froissart disk 
. But there, with the growth of energy, more higher transverse momenta scales enter the 
game. And although on every such scale saturated density is finite, the total density grows 
with E and the disk becomes more and more black. How this question is solved for gravity is 
unclear. 
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Usually in this case one starts from a Green function corresponding to a ba- 
sic Regge object (like bare pomeron) and then, using it, one can build various 
higher order diagrams for amplitudes and inclusive cross-sections. In our case 
this implies to start from 2G pole with a(0) = 3 + A corresponding to (p3]). But 
we make one step forward and start directly from the "economized" 2G object. 
For pomerons this step corresponds to a transition to "Froissarons" , which are 
the economized sums of pomerons, and then to construction of higher diagrams 



using such objects 24 . Applying this approach to our case we form the basic 
reggeon Green function - proportional to the amplitude corresponding to the 
S-matrix @ : 

D(y,b) = i(l-S(y,b)) ~ i 9{R 2 ± (y) - b 2 ) + D s °f\y,b) , (36) 

where D so; ' t (y,b) is the soft part of D containing gravitons with ui, k± <C 
m p . Using this D(y, b) one can construct higher reggeon diagrams and various 
inclusive cross-sections. The hard inelastic cross section is given by 

a m (Y) = J d 2 B lmD(Y 1 B) ~ nR\{Y) + <r£ ft (Y) (37) 

as considered in the previous section. The corresponding inclusive cross-section 
for hard gravitons with the definite impact parameter b and the rapidity y and 
the initial Y and B is given by the reggeon diagram with two cut D s - it takes 
in the impact parameter representation a simple form 

p(y,b,Y,B) ~ V1mD(y,b)-ImD(Y-y,B-b) = (38) 
= V8(Rl(y)~b 2 )-9(Rl(Y-y)-(B~b) 2 ) + p so ^ t (y,b,Y, B) , 

where V - is the inclusive graviton emission vertex, and the p^f 1 - the soft con- 
tribution. Other higher inclusive cross sections can be constructed in the similar 
way. The region where such p(y, b, Y, B) is big is shown in Fig. 4. The higher 




Fig. 4 : The shape of region of intersection of black disks in 
(x± * y) space at various B. The width of shaded region at 
given rapidity is connected with value of inclusive cross-section as 
represented in Fig. 6 



loop corrections in D do not change the structure of the inclusive cross-section 
in the "black" part of D(y,b) i.e. at b < R±(y), but can be essential at the 
border of the AS-disk, where their influence on a diffractive and a soft radiation 
processes (like in the case of Froissaron p3 ) is not a small. 
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Now let us try to interpret the expressions ( |3S| ) in the parton language. In 
terms of multiparticle Fock space components this is rather a complicated task 
and is not yet considered in full details even for a simple scalar multiperiphcral 
interaction. But we know many qualitative aspects of such processes from the 
space-time interpretation of various ladder and eikonal diagrams corresponding 
to a reggeon exchange. We consider separately soft and hard gravitons. Because 
most of partons in the AS disk have energies uj ~ m p we will consider firstly 
only them. 

Hard graviton production. 

At given rapidity y and x±_ hard gravitons can be produced only from a collision 
at the same rapidity and x± of black components of two disks. This corresponds 
directly to the "two 6" intersection term in (^8|). Then the full inclusive spec- 
trum at the same y is proportional to the area of the intersection of two black 
parts of AS disks at the given rapidity. So, for definite B and y, we can have 
three different configuration of disks intersections as shown in Fig. 5. In the 




Fig. 5 : Configuration of intersecting black AS disks is sown at 
three rapidities: A\ - closer to the fragmentation region, A3 - in 
the center of mass system. This can be compared with spectra at 
Fig.6 at B > B cr 

cases (^4i) and (^2) the two AS-disks intersect - and the area of the intersec- 
tion, in Planck units, A(B, R±(Ei), R±(E 2 )) defines the value of hard inclusive 
cross-section. In the configuration (b) there is no intersection of the black parts 
of the disk - so the contribution to the hard inclusive cross-section at these B 
and y is zero. We neglect effects of smeared hard disk borders and then the 
hard inclusive density at given B and y is approximately given by 

p{B,y,Y) ~ m 2 p A(B,R u R 2 ), (39) 

where R\ = R±(Ei), R\ = R±(E\), E\ = m p e v , E2 = m p e Y ~ v , and the area 
of the black disk intersection can be represented as 

A~ J d 2 Xl _ 6{R\-x\) 6{Rl- \x±~B\ 2 ) (40) 

A = CiRf + C2RI for \Ri - R 2 \< B < R l + R 2 , 

A = min(7ri?f, ttR 2 2 ) for B <\R X - R 2 \ , 

where 

Ci = 0i - 5 sin2 #i . C2 = O2 - \ sin20 2 , 
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2BR t ' 2BR 2 

The structure of hard graviton spectra in rapidity, corresponding to (^ , ^p| ) , is 
illustrated at Fig. 6 It changes from highly peaked in c.m.s. for B ~ to 
the purely diffractive like spectra with two filled bumps near the fragmentation 
region of colliding particles for large B ~ m~ 2 \fs. The gap in inclusive spectra 




Y/2 V y Y/2 



Fig. 6 : Hard gravitons inclusive spectra in rapidity (a,b,c,d) 
for various values of the impact parameter B, where B cr = 

around the central position y = Y/2 is formed when B > B cr , where the 
B cr correspond to a configuration when two hard disk in c.m.s. only touch 
each other, i.e. 2R±(Y/2) = B cr . The multiplicity of produced gravitons and 
the maximal value of the spectra at Fig. 6 also depends strongly from B, and 
drastically changes their behaviour at B ~ B cr ~ m" 1 (s / rn^) 1 / A . From ( ftO[ ) 
one can simply estimate the behaviour of the mean multiplicity of hard gravitons 
as a function of s and B: 

ma,.) - .(«,-*) g/G^j) + •(»-*-) sjp. <«) 

where the function f(B/B cr ) = m v s~ x l' 1 j dyA(B, Rj_(y), Rx(X — y)) depends 
slowly on the argument. The multiplicity distribution due to ( |40| ) is also mainly 
defined by geometrical parameters of collision, and the probability to produce 
n hard gravitons 

1 /to 2 \ -3/2 m 

w n ~ 4 + (— ) h{n^p)Q{n - n cr ) (42) 

71 c> \ S J y/S 

is fast decreasing almost for all n. The small additional term with the slowly 
varying function f\ ~ 1 in r.h.s. of ( ^2[ ) is essential only for n > n cr ~ ^fsjra v . It 
comes from collisions with B < B cr (s). Although the total inelastic cross-section 
( [37| ) for the production of hard gravitons is big ~ sm~ 4 , the corresponding 

12 Later on we put A = 0. This simplifies expressions and at the same time corresponds 
probably to the correct value. 
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final state contains a small number of particles (Fig. 6 d). The processes (Fig. 6 
a,b,c) , where the hard multiplicity in final state is high ~ t/s, come only from 
collisions with B < B cr - their cross-section is ~ m p 3 ^ 2 y / s. 

The mean hard graviton multiplicity corresponding to (^2|) is low < n >~ 1 
at all s, and the dispersion of the multiplicity grows: < n— < n >> 2 ^ Ins, but 
not so fast 

Note that as a result of such behaviour of w n the full inclusive spectra of 
hard gravitons integrated over B takes a very simple scaling form 

P(V, Y) = err 1 (Y) ( d 2 B p(B, y, Y) ~ (43) 



™% Rj R\ E 1 ■ E 2 
- ^271 = n 5 const(y,Y)~l 

like in the case of constant asymptotic cross-sections. Although here the fluctu- 
ations in shape of the individual events can be very big and are mainly governed 
by the distribution of impact parameters, like in the case of a nucleus-nucleus 
collision. 



The role of soft gravitons and the longitudinal boost behaviour. 

We have not considered explicitly the production of low energy gravitons with 
u> <C rn p and of the soft gravitons with k± <C m p . There are two types of 
such processes. In one, which takes place on large impact parameters, only the 
soft gravitons are produced. Such processes must be more or less satisfactorily 
described by the perturbative methods - and the complete space-time picture 
will be not far from that corresponding to a classically described gravitational 
wave radiation during a fast particles collision, for example like that given in Q . 

The other processes, taking place at smaller B, also contains soft gravitons 
accompanying hard gravitons created as a result of a collision of black AS disks. 
But here the soft and low energy partons also enter directly in the production of 
hard gravitons. The consideration of their interaction is needed to consistently 
explain (in parton terms) the generation of hard particles for relatively big 
impact parameters B cr < B < m p 2 ^fs , and the inclusive spectra of type 
Fig.6d, when we view on a collision from the c.m.s., where the hard disks do 
not collide directly. 

The related more general question concerns the invariance of cross-sections 
and reactions outcomes under longitudinal boost transformations. The parton 
state of fast particles changes under such boosts in a complicated and nontriv- 
ial way: the mean number of partons and their space configuration changes 
etc. Therefore the requirement of the invariance of various cross-sections calcu- 
lated in the parton approach under boost, corresponds to very strong condition, 
essentially restricting the structure of the parton state itself FJ. 

13 This radically differs from possible expectations, when for all B < nip 2 ^/s the big black 
hole (possibly rotating) is created, leading at the end to a almost full dissipation of the initial 
energy into soft in c.m.s. particles. 

14 In the massive theories like QCD, the rate of the growth of a parton black disk ra- 
dius R±(E) is almost completely fixed by the invariance of the inelastic cross-section 
ffj„ = tt(R±(Ei^) + R±(E2/£,)) 2 under changes of the boost parameter §. Such a condi- 
tion leads only to two solutions : Rj_(E) = const(E) and R±(E) = c\nE. This corresponds 
either to ai n (E) = const or ai n (E) ~ In 2 E - that is the Froissart-like behaviour. Moreover, 
by same method one can show that in the Froissart case only the black (not a gray) disk case 
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Let's consider the interaction of the AS disk with E\ ^S> m p with an almost 
massless particle which can have a "low" energy E% <§; m p - so that it can be 
considered as "test" objects and in the same way choose the full s ~ 2E\E 2 ^> 
m 2 p . 

Then examine the behavior of this system under the longitudinal boost - 
when the transformed energies become 

E x - E 1 $ , E 2 -> E 2 H , (44) 

but the impact parameters B and s are not changed. 

Firstly choose B = B\ < R±(s) ~ ™p 1 ( s / m p) 1 ' 4 so that the hard gravitons 
are produced in the final state. If £ are such that E 2 ~ m p then we have the 
collision of the ASi black disk with a hard particle, which is inelastically scat- 
tered by a disk with the probability = 1. After that the production of secondary 
hard gravitons starts almost immediately and the real particles appear at times 
~ m~ . Call this frame - frame I and fix £ ~ 1 in this frame. 

If we choose next £ such that E 2 -C m p then we must suppose that the 
interaction of the U E 2 " particle with the hard saturated AS\ disk at the same 
B takes place with the same probability = 1, although the cross-section of 
the interaction with individual partons decreased ~ l/£ 2 - Now the AS\ parton 
system is only "softly" excited during the interaction. After that the instabilities 
in the AS\ disk grow gradually and first hard gravitons will be created only after 
the AS disk moves in z on ~ m" 1 £ from the collision place. Let's call this frame 
- frame II. 

But for such £ (in frame II) the radius of the black part of the ASi disk, 
where partons are in the saturated phase, is also increased 

JZlCEi) -» R±(Ekt) = ti 1/2 R±(Ei) . 

Therefore the "E 2 " particle should be absorbed by a disk with the probability 
= 1 for larger B up to ^ 1 ^ 2 R±(Ei). This probability of the particle U E 2 " interact 
at B — B 2 where R±(E\) < B 2 < R±(Ei£) must not depend on £. In the frame 
I at B = B 2 the 11 E 2 " particle does not collide directly with a black disk. Now 
directly with the black disk can interact only soft partons accompanying the 
particle U E 2 " . And this mechanism can regulate the boost invariance of the 
probability interaction = 1. This role of soft partons-gravitons is even more 
evident if, at the same B = B 2 , we select £ so to move to cm. system, where 
the black disks are most far one from another in the transverse direction. 

One can estimate simply maximal impact parameters for which the capture 
of soft partons from one colliding particle by the black disk of another particle 
take place with the probability = 1. For the capture of soft partons we need to 
choose B and u) so that as minimum as one soft parton from ll AS 2 " falls into 
the black disk of ASi or vice versa - this gives : 

n(E 2 ,uj, B) ■ R 2 l (Ei 1 lu 1 ~ m p ) ~ m~ 2 

From the other hand, the size of the transverse localization of the soft parton, 
which is ~ must be much less than the AS\ black disk radius R± {E\ , u)\ ~ 

is allowed. The other restriction can also be found by the same way. The reason why the bust 
invariance of cross-sections is so restrictive for parton states is probably connected with that 
it plays the role of t-unitarity which is very essential for high energy amplitudes, but can not 
be explicitly imposed on parton wave functions W\- 
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m p ). Otherwise the soft parton will mainly elastically scatter on the black disk 
and not be captured by him. This corresponds to the condition for minimal u) : 

lu ■ R±(Ei, m p ) ~ 1 

Then, using expressions ( |33| ) , ( p5| ) for n and R± , it is simple to find B — Rj_ ma x 
and uj fulfilling these two conditions in the case of the arbitrary intercept A : 



(2+A)/4 /BAHfA+^A 2 /£ 2 \l+A/2 



m 



The inelastic processes up to such big impact parameters B ~ R±max are es- 
sentially influenced by "strong" gravity mechanisms. The corresponding hard 
inelastic cross-section 

„ / R, \ 1+fA+iA 2 , R s l+A/2 

<y m {E u E 2 ) ~ (i? ±mQ:c ) 2 ~ m p 2 (±1) f^ 2 -) (45) 



can be boost invariant under (44) only if the powers of E\ and i?2 in (|45|) are 
equal - because only in this case 2E\E2 ~ s does not depend on £. This gives 
condition on A with two solutions 

A = , A = - 2 

The first solution corresponds to the hard <Tj„ ~ sm~ 4 cross-section which we 
have discussed above and gives the additional arguments for the A = case. 
The second solution corresponds (a little peculiarly) to asymptotically constant 
cross-sections, typical for hadron interactions and associated with a pomeron. 



5 Instability of the final state. 

The black holes creation after collision 

Before the collision, in parton black disks of incoming particles, we already 
have too much 'low' energy transverse partons-gravitons so that in fact the 
full virtual energy is concentrated in them. But this state is evidently stable 
under the implosive collapse due to strong phase correlations between these 
partons. During the collision some partons are exited and "discorrelated" , and 



later become free with the inclusive spectra of type (38). Only after that the 
possible long-range gravitational instabilities in the evolved particles system can 
lead to the formation of various density clusters over the "initially smooth" p(y) 
background. 

The main attractive instability can take place for groups of created parti- 
cles having close longitudinal velocities. We can divide conditionally the p(y) 
spectrum in layers in the rapidity of the width i5 ~ 1 in such a way that the rel- 
ative energies of particles(gravitons) from the same layer are not high, and the 
hard partons from the neighbor layers look ultrarelativistic. In a first moment 
the created particles density in every such layer is high - close to the critical. 
Therefore these particles (gravitons) strongly interact one with another during 
some period of time before there density becomes so small that they can move 
almost free. During this period of time the relative energies of particles in layer 
are partially thermalized and they behave in fact as a massive particles. It is 
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essential that because of such multipole scattering on "media" , for some period 
of time particles can not escape from layer as free massless particles. In fact such 
condition fixes the value of S and the portion of particles that are "confined" in 
layer. 

If we go to a longitudinal system with a definite rapidity yi, where some 
layer is "standing" (marked in Fig. 7) then we see that particles from other 
layers move fast away. In fact in the corresponding time scale these particles 
can be still in the virtual state and become free only after time t ~ exp \y — y±\. 
So we can consider this layer separately from the others. Its mass is 



P(y) 




Fig. 7 : Inclusive spectra for hard gravitons divided on layers in 
rapidity. These layers can transform during the following evolu- 
tion into black holes with the same rapidities and masses ~ p(y). 

M( yi ) ~ m p A(y 1: Y,B) ~ p{ Vl ,Y,B) , (46) 

where A is the transverse area of the "yi" layer. This mass is concentrated in the 
region of the space with the longitudinal size ~ m^ 1 exp (8) and the transverse 

size ~ y/A. These magnitudes are much less than the size of the horizon of 
the BH with the same mass M(yi), and at the same time particles constituting 
this mass are not ultrarelativistic. Therefore the most part of particles from 
such a layer at time ~ \f~A~ transforms in a black hole with the mass M(y±) 
and the transverse position, defined by the center of the region A. During this 
time the horizon "surrounding" these particles transforms from a pancake-like 
configuration to sphere. 

The neighbours to "t/i" layers u y\ — 5" and a y\ + 8" also transform to black 
holes with masses m p p(yi — 6) and m p p(yi + 6) , but this will take place later 
and they will have the relative to "j/i" velocities ~ 1. So gradually the chain 
of black holes will be produced whose distribution in longitudinal momenta is 
typically multiperipheral and close to the uniform p*) - with p = const(y) ~ <5 _1 . 
But their masses are different at fixed B and are given by values of p(yi, Y, B). 
On average - when we integrate over B then the mean masses of produced BH 
are also constant over y and small ~ m p , because the main contribution comes 
from the large B. 

Because every such BH is produced from a layer of partons having approx- 
imately the same longitudinal rapidity, these BH will be non-rotating. This is 
accurate within the small "thermal" fluctuation, and border effects. Because 
the borders of layers are slightly asymmetric for B ^ (see Fig. 4) some angular 

15 In fact uniform distribution in rapidity can be slightly distorted, because the effective 
widths of layers in rapidity can be correlated with there masses. 
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momentum can be concentrated there. The main part of the angular momen- 
tum from the initial state (especially for large B) is transferred to the relative 
motion of the created BH and to the surrounding particles not captured by BH. 

The rapidity distribution of produced BH in the individual events is con- 
trolled mainly by the corresponding values of B and repeats the curves in 
Fig. 6. For B < B c = m~ x (E/rrip) 1 ^ we have a continuous (multiperipheral) 
in y chain of BH, whose masses are given by curves Fig.6(a-b). For B > B c two 
chains of black holes moving in forward and backward directions are created - 
their masses repeat again the behaviour of p(y, Y, B) in Fig.6(c,d). These chains 
are separated by the rapidity gap of width Ay ~ Y — 2\n(Bm p ) and for big 
Ay when B ~ y/s/rrip correspond to a typical diffractive generation of black 
holes [3 

Probably the part of particles falling between layers can escape and not 
captured by one of BH. Their percentage is possibly not big, but it is complicated 
to estimate them by the qualitative methods used here. 

It remains a question : can some of these black holes (or all of them) merge 
during the future evolution, if the rate of longitudinal grow of there horizons 
is faster than these BS separate one from another. In fact this is a purely 
classical GR problem, and it needs additional investigation. Here we present 
only following qualitative argument. 

The longitudinal structure of particles system created in a collision of black 
AS disks is similar to the one-dimensional homogeneous expanding cosmological 
solution. In both cases particles are "injected" initially in such a way that in 
every longitudinally boosted frame we locally see particles flying away with the 
same 2D density and the momentum distribution. Usually for such expanding 
system the small long-range density perturbations are damped with time or not 
growth. Only a short-range (scale S) large density fluctuations, created on the 
initial "viscous" stage, can growth and clusterize. But on a scale much larger 
than this no new structures emerge. 



6 Conclusion and the final remarks 

Thus we come to a picture of gravitational interactions at E 3> m p which is 
in some respect reminiscent to the picture of asymptotic hadronic (QCD) in- 
teractions in the Froissart limit: the transplanckian particles collision looks like 
an interaction of black disks with the production of particles in the geometrical 
intersection region of disks. The essential difference is that for the gravitational 

16 Above we discussed how complicated can be the "trajectories" leading from the colliding 
particles state to a final state containing the BH. Such trajectories are evidently absent in the 
semiclassical probability estimate exp ( — S gr (g)^j , when we put in Euclidean classical action 
S gr some continued solutions g containing the BH. There are two ways to cure this. In one, we 
must put all higher radiative correction in the Euclidean action (the necessary number grows 
with the energy). Then the corresponding equation of motion can have necessary trajectories 
g with a low value of effective action and a high statistical weight. The other way is to use 
the classical action and consider the transition to the BH not from the few particle initial 
state, but from the specially prepared multiparticle state (parton component), containing 
approximately the same number of degrees of freedom as the entropy of the BH. But in both 
cases the "advantages" of the Euclidean approach disappear and the consideration in the 
Minkowski frame is much simpler. This is in fact the general reason because the Euclidean 
methods do not work properly at high energies. 
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interaction the radiuses of these black disks grow much faster ~ E 1 / 2 with en- 
ergy, as compared to the QCD case, where these radiuses grow relatively slow 
~ In E. The other distinction is the character of the interaction in the final 
state. In the QCD the final state hadronization is short-range and so this does 
not essentially change the structure of the final state. But for the gravitational 
case we have a long-range attractive interaction which can create a "multipe- 
ripheral" chain of black holes. Finally all this difference is the consequence of 
the masslesness of gravitons. If QCD were massless at large distances then prob- 
ably also there can some similar phenomena originate - like powerlike growth 
with energy of inelastic cross-sections, etc. 

Additional remarks in conclusion. 

• At "not too high" energies s ~ m 2 , when both black disks from hard 
partons are only in the embryonic state, it is complicated to make any 
quantitative conclusion (in particular - predict the threshold for the first 
BH creation and their mass) by the methods we used. Here all essential 
parameters are of order of unity and the possible dynamics can be rather 
involved. But there are also no reasons to expect that the purely classical 
estimates give reliable values. 

• This model can probably extended to the case with additional hidden 
dimensions D > 4. Here one can distinguish two cases 

a) "Small" additional compact dimensions, with the radius close to the 
Planck scale ~ m p . In this case all higher Kalusa-Klein modes act simply 
as massive particles strongly coupled at the scale ~ m p to gravitons. Prob- 
ably they do not change any qualitative aspect of picture, but contribute 
"only" to various renormalization of numerical parameters. 

b) If the size of additional dimensions (or of some of them) in which 
gravitons can penetrate is very large compared to the 4 D Planck scale, as 
in the case of the TeV gravity O], then some aspects of the picture can 
change more essentially. In this case too much depends on the concrete 
realization of the gravity in the brane world. For example, in the case of 
only large and non warped additional dimensions we simply have the D > 
4dim. gravity, but some particles (and the colliding ones also) are glued 
to 3 space dimensions. Then we have the transverse graviton spectrum 
dn ~ (E 2 duj/uj 3 )d D ~ 2 (k±/mp). This leads to a black D-2 dimensional 
disk filled with partons at the Planck density with the radius 

/ \ l/(£>— 2) 

R(E) ~ m-UE/mpJ 

and to cross-sections <7j n ~ R 2 (E). By the same way one can find the 
created hard graviton density p(y,B,Y ~ mp~ 2 A((B,y,Y)), where A 
is now the D — 2 dim. volume of the colliding disk intersection region. 
Qualitatively the B dependence on p is the same as in Fig. 6. In the final 
state one can again expect the gravitational long-range instability and the 
produced particles clusterization into the chain of black holes. But now 
(and here much depends on the details of brane world models) created 
gravitons and black holes can move in large additional dimensions and 
dissipate there. In corresponding "accelerator events" this can look out 
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as if a part of the colliding energy escapes. Evidently there can be many 
variations of this simple scenario. 

• In this model the elastic scattering amplitude contains components coming 
from different sources. One is Coulomb like - it comes from the "classical" 
soft graviton exchange at large impact parameters. Of same type are the 
soft long-range multiperipheral corrections. The other one corresponds to 
the imaginary diffraction contribution to the elastic amplitude generated 
by the big inelastic cross-section ~ m~ 2 {s / m^) . 

The finite angle high energy elastic scattering usually comes from small 
impact parameters and is triggered by the fluctuations in which in an 
initial state one has only few (minimally possible) number of partons, 
concentrated in a small size volume. 'For transplanckian collisions this 
corresponds at least to components of partonic wave functions of both 
particles without hard black disks. This probability is ~ exp(— s/rrip C), 
and this leads to the too fast decreasing term in the finite angle scattering 
cross-section : 

da/dO S -j expf — %C(6)) + (da/dO) 

nip V m£ J \ J soft 

where \da/dO\ ~ exp(— ciy/s ) is the large impact parameter con- 

V / soft 

tribution, generated by the multiple scattering. This behaviour can be 
compared with the Cerelus-Martin bound and with the finite angle cross- 
section for a string scattering p5| . 
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